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Abstract 

Reliable  predictions  for  hypersonic  laminar- 
turbulent  transition  will  have  to  be  based  on 
simulation  of  the  transition  mechanisms.  Accurate 
simulation  of  the  second-mode  instability- wave 
growth  on  a  round  cone  at  zero  angle  of  attack 
would  be  a  significant  step  toward  this  goal.  Hy¬ 
personic  stability  experiments  in  this  geometry 
are  reviewed,  and  three  are  selected.  Extensive 
and  previously  unpublished  data  from  Stetson’s 
sharp  and  blunt  cone  experiments  are  reported. 
Approximate  computations  were  carried  out  using 
inviscid-boundary-layer  and  eN  methods.  The 
experimental  results  are  compared  to  previous 
high-accuracy  computations  and  to  the  present 
approximate  computations.  Quantitative  agree¬ 
ment  is  at  present  insufficient  for  reliable  eN 
predictions.  The  difficulties  with  both  experiment 
and  computation  are  described. 

Introduction 

General  Issues 

Laminar-turbulent  transition  in  high-speed 
boundary  layers  is  important  for  prediction  and  con¬ 
trol  of  heat  transfer,  skin  friction,  and  other  bound¬ 
ary  layer  properties.  However,  the  mechanisms  lead¬ 
ing  to  transition  are  still  poorly  understood,  even  in 
low-noise  environments  [1].  Applications  hindered 
by  this  lack  of  understanding  include  reusable  launch 
vehicles,  high-speed  interceptor  missiles,  hypersonic 
cruise  vehicles,  reentry  vehicles  (RV’s)  [2],  and  dis¬ 
crimination  among  incoming  RV’s  [3]. 
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The  transition  process  is  initiated  through  the 
growth  and  development  of  disturbances  originating 
on  the  body  or  in  the  freestream  [4].  The  recep¬ 
tivity  mechanisms  by  which  the  disturbances  enter 
a  boundary  layer  are  influenced  by  roughness,  wavi¬ 
ness,  bluntness,  curvature,  Mach  number,  and  so  on. 
The  growth  of  the  disturbances  is  determined  by  the 
instabilities  of  the  boundary  layer.  These  instabili¬ 
ties  are  in  turn  affected  by  all  the  factors  determining 
the  mean  boundary  layer  flow,  including  Mach  num¬ 
ber,  transverse  and  streamwise  curvature,  pressure 
gradient,  and  temperature  [1].  Relevant  instabilities 
include  the  concave- wall  Gortler  instability  [5],  the 
first  and  second  mode  TS-like  instability  waves  de¬ 
scribed  by  Mack  [6] ,  and  the  3D  crossflow  instability 
[7].  The  first  appearance  of  turbulence  is  associated 
with  the  breakdown  of  the  instability  waves,  which 
is  determined  by  various  secondary  instabilities  [8]. 
Local  spots  of  turbulence  grow  downstream  through 
an  intermittently-turbulent  region  whose  length  is 
dependent  on  the  local  flow  conditions  and  on  the 
rate  at  which  spots  are  generated  [9]. 

In  view  of  the  dozens  of  parameters  influencing 
transition,  classical  attempts  to  correlate  the  general 
transition  ‘point’  with  one  or  two  parameters  such  as 
Reynolds  number  and  Mach  number  can  only  work 
for  cases  that  are  similar  to  those  previously  tested. 
For  general  flight  data,  the  scatter  is  large  (e.g.,  Ref. 
[10]).  Transition-estimation  methods  that  are  reli¬ 
able  for  a  broad  range  of  conditions  will  need  to  be 
based  on  an  understanding  of  the  physical  mecha¬ 
nisms  involved. 

The  simplest  and  best  developed  of  the 
mechanism-based  methods  are  the  eN  methods, 
which  attempt  to  correlate  transition  with  the  in¬ 
tegrated  growth  of  the  linear  instability  waves.  Al¬ 
though  these  methods  neglect  receptivity  and  all 
nonlinear  effects,  such  as  wave  interactions,  nonlin¬ 
ear  breakdown  effects,  roughness,  and  so  on,  they 
have  shown  promising  agreement  with  experiment 
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[1].  Agreement  is  fairly  good  for  a  variety  of  con¬ 
ditions  where  the  environmental  noise  is  generally 
low  [11].  However,  the  accuracy  of  many  high-speed 
eN  computations  remains  uncertain,  for  it  is  often 
not  clear  whether  bluntness  effects,  wall  temper¬ 
ature  distributions,  and  so  on  were  handled  with 
sufficient  accuracy  [12].  Highly  accurate  mean-flow 
profiles  are  essential,  since  the  stability  equations 
are  very  sensitive  to  the  second  derivatives  of  these 
profiles.  Although  wave  interaction  effects  and  3D 
effects  can  sometimes  be  handled  with  correlations 
[13],  eN  methods  are  only  an  intermediate  step  along 
the  way  to  reliable  mechanism-based  methods. 

Direct  simulations  of  transition  [14]  and  the 
recently  developed  Parabolized  Stability  Equations 
(PSE)  [15]  have  in  some  ways  advanced  theoretical- 
numerical  work  far  ahead  of  the  experimental 
database.  The  numerical  work  is  not  yet  able  to 
include  complex  effects  such  as  roughness,  wavi¬ 
ness,  internal  shocks,  and  most  bypasses.  Three- 
dimensional  mean  flows  and  their  instabilities  are 
only  beginning  to  be  treated  correctly  from  first 
principles  [16].  The  simulation  of  bluntness  effects 
[17]  and  chemistry  effects  [18,  19]  from  first  princi¬ 
ples  is  only  beginning.  When  the  numerics  are  based 
on  the  correct  physical  mechanisms,  however,  they 
can  provide  much  more  detail  regarding  the  transi¬ 
tion  process.  Experimental  work  that  describes  not 
only  the  location  of  transition  but  also  the  mecha¬ 
nisms  involved  is  needed  in  order  to  improve  these 
modern  theories.  The  key  mechanisms  need  to  be 
identified,  in  part  through  experimental  work,  and 
the  key  numerical  results  need  to  be  validated  ex¬ 
perimentally. 

Unfortunately,  most  of  the  ground-test  data  are 
ambiguous,  due  to  operation  in  high-noise  conven¬ 
tional  wind  tunnels  and  shock  tunnels,  with  dis¬ 
turbance  levels  much  higher  than  in  flight  [20]. 
The  mechanisms  of  transition  operational  in  small- 
disturbance  environments  can  be  changed  or  by¬ 
passed  altogether  in  high-noise  environments  [21]. 

Transition  on  Reentry  Vehicles 

Laminar-turbulent  transition  is  the  dominant 
error  source  in  aerothermal  heating  analyses  of  bal¬ 
listic  reentry  vehicles.  This  is  shown,  for  example,  in 
Fig.  8  of  Ref.  [22] ,  which  presents  computations  and 
measurements  of  the  surface  heat  transfer  during  the 
Reentry-F  test.  Agreement  is  good  for  both  the  lam¬ 
inar  and  turbulent  regions,  once  the  transition  loca¬ 
tion  is  known.  Harris  Hamilton  from  NASA  Langley 
says  that  typical  accuracies  are  20-25%  for  the  tur¬ 
bulent  boundary  layer,  and  15-20%  for  the  laminar 


layer.  The  largest  uncertainty  is  due  to  transition, 
which  causes  a  large  increase  in  the  heat  transfer 
rates  (private  communication,  March  1999).  Similar 
results  are  reported  in  Ref.  [23].  Ref.  [10]  shows  typ¬ 
ical  empirical  correlations  such  as  Reg/Me ,  around 
which  the  data  scatter  by  as  much  as  a  factor  of  3. 

Real  ballistic  RV’s  are  ablating,  rough,  and  at 
non-zero  angle  of  attack.  Ground-test  measurements 
of  transition  will  not  be  able  to  simultaneously  simu¬ 
late  all  the  characteristics  of  the  high-enthalpy  flight 
environment;  thus,  ground-based  experiments  must 
discover  and  document  the  transition  mechanisms, 
in  order  to  develop  and  validate  computational  mod¬ 
els  which  are  then  extended  to  flight  conditions 
[24].  For  ground-test  measurements  of  stability,  a 
smooth-wall  perfect-gas  cold-flow  round  cone  with 
small  bluntness  is  a  reasonable  first  approximation. 
Measurements  of  the  instability  waves  (or  transition 
mechanisms)  are  needed  under  these  conditions,  for 
comparison  to  computation.  A  number  of  stabil¬ 
ity  experiments  have  been  carried  out  for  this  case. 
Although  all  contain  significant  flaws,  a  detailed  re¬ 
analysis  of  this  data  is  being  carried  out  for  NATO 
RTO  Working  Group  10  (WG  10).  The  present  work 
summarizes  this  reanalysis,  and  presents  compar¬ 
isons  to  computations  reported  in  the  literature,  and 
to  limited-accuracy  eN  computations  carried  out  by 
the  author. 

Review  of  Stability  Experiments 


A  necessary  step  in  the  further  development  of 
mechanism-based  methods  is  accurate  computation 
of  the  actual  instability-wave  growth.  At  high 
speeds,  this  is  a  major  challenge.  At  low  speeds, 
such  computations  were  slow  to  become  validated 
(e.g.,  Ref.  [25]);  this  validation  is  more  critical  at 
high  speeds,  where  there  are  many  more  physical 
factors  whose  influence  must  be  properly  taken  into 
account.  However,  even  in  noisy  conventional  tun¬ 
nels,  there  are  few  measurements  of  the  instability 
mechanisms  causing  transition.  Of  these  instabil¬ 
ity  measurements,  few  have  been  carried  out  with 
calibrated  instruments  and  well-documented  condi¬ 
tions,  so  that  a  fairly  reliable  comparison  to  com¬ 
putations  can  be  carried  out.  Rather,  existing  sta¬ 
bility  experiments  have  mainly  served  to  discover 
the  stability  phenomena,  confirm  the  fundamental 
aspects  of  the  theoretical  predictions,  and  supply 
preliminary  quantitative  information.  While  such 
experiments  have  been  an  important  step  forward, 
the  requirements  for  code  validation  experiments  are 
much  more  stringent  (e.g.  Refs.  [26]  and  [27]).  Even 
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for  the  surface  pressure  distribution  under  cold-flow 
conditions,  accurate  validated  results  remain  to  be 
obtained  [28]. 

The  transition  members  of  WG  10  felt  that 
the  following  elements  should  be  present  in  code¬ 
validation  experiments  for  hypersonic  transition: 

1.  Detailed  and  reliable  measurements  of  the  tran¬ 
sition  mechanisms. 

2.  Accurate  knowledge  of  the  mean  flow.  Re¬ 
peated  checks  of  boundary  layer  symmetry,  in¬ 
strumentation  calibrations,  tunnel  flow  nonuni¬ 
formity  effects,  repeatability,  etc.  Accurate 
agreement  with  computation  of  mean  flow.  This 
requires  repeated  tunnel  entries  and  coopera¬ 
tion  between  experiments  and  computations. 

3.  Linear  instability  comparisons  require  cali¬ 
brated  measurements  of  the  fluctuations  at  a 
known  position  in  the  eigenfunction.  It  seems 
possible  that  these  might  be  carried  out  in  con¬ 
ventional  tunnels  using  ensemble  averaging  of 
controlled  disturbances. 

4.  Nonlinear  secondary  breakdown  effects  are  de¬ 
pendent  on  small  fluctuations  combining  with 
the  primary  instability.  It  seems  doubtful  that 
these  can  be  repeatably  and  successfully  stud¬ 
ied  except  in  quiet  tunnels,  since  even  controlled 
secondary  disturbances  may  be  swamped  by 
tunnel  noise.  For  similar  reasons,  receptivity 
experiments  may  also  require  quiet  flow. 

These  specifications  are  updated  from  Ref.  [29], 
which  gives  additional  suggestions,  and  shows  the 
long-term  importance  of  a  coordinated  approach. 

Unfortunately,  these  specifications  are  ex¬ 
tremely  difficult  to  meet  at  hypersonic  speeds.  Al¬ 
though  new  experimental  efforts  are  underway  at 
Purdue  and  at  ITAM,  Novosibirsk,  Russia,  all  ex¬ 
isting  data  fall  short.  Despite  these  shortcomings,  it 
seemed  valuable  to  assess  the  best  existing  datasets 
in  detail,  as  both  a  determination  of  the  state  of 
the  art,  and  as  a  means  of  finding  the  improvements 
that  are  needed.  Three  datasets  were  identified  for 
further  analysis.  New  experimental  measurements 
are  presently  being  obtained  at  ITAM,  Novosibirsk, 
Russia  [30];  these  are  not  complete,  and  have  not  yet 
been  examined  in  detail. 

While  the  following  three  datasets  all  contain 
significant  flaws  and  limitations,  they  are  the  best 
presently  available,  in  the  opinion  of  the  team.  Ac¬ 
curate  computation  of  the  wave  growth  and  eigen¬ 
functions  for  these  cases  would  be  a  significant 


step  towards  the  further  development  of  mechanism- 
based  methods.  Comparisons  between  the  experi¬ 
ments  and  the  results  of  multiple  independent  com¬ 
putations  would  shed  light  on  the  strong  and  weak 
elements  in  all,  and  aid  in  planning  new  work. 

Because  of  the  limited  resources  presently  avail¬ 
able,  attention  is  focused  on  axisymmetric  geome¬ 
tries.  The  computational  resources  required  are 
much  less  than  in  the  3D  case.  All  of  the  datasets 
are  for  cold  hypersonic  flow  without  chemistry,  be¬ 
cause  detailed  measurements  of  the  mechanisms  of 
high-enthalpy  hypersonic  transition  do  not  exist.  Al¬ 
though  the  data  shown  below  are  primarily  measure¬ 
ments  of  the  cold-flow  second-mode  instability,  other 
mechanisms  do  of  course  remain  important. 


Present  Computations  Using 

Approximate  Methods 

Several  of  the  experimental  flows  were  computed  by 
the  author  using  fairly  elementary  methods  that  are 
commonly  available  to  nonspecialists.  The  inviscid 
axisymmetric  flows  were  computed  by  Sandia  Na¬ 
tional  Labs  and  provided  by  Dr.  Dave  Kuntz.  The 
2IT  code  was  used  to  compute  the  inviscid  nosetip 
flow  [31].  This  code  solves  the  unsteady  Euler  equa¬ 
tions  using  a  time-dependent,  finite-difference  tech¬ 
nique,  with  the  cylindrical  flowfield  mapped  into 
a  computational  domain.  The  technique  is  shock¬ 
fitting  in  nature,  with  the  characteristic  compatibil¬ 
ity  relations  invoked  at  the  bow  shock  and  imper¬ 
meable  surface  boundaries.  MacCormack’s  explicit 
scheme  is  used  to  integrate  the  governing  equations, 
producing  a  solution  technique  which  is  second  order 
in  both  time  and  space.  The  solution  is  converged 
in  time  to  produce  both  a  steady-state  flowfield  so¬ 
lution  on  the  spherical  nosetip  and  an  initial  data 
plane  for  the  afterbody  code. 

The  afterbody  flowfield,  surface  pressure,  and 
shock  shape,  were  solved  with  the  SANDIAC  code 
[32],  an  extension  of  the  GE-3IS-SCM/ACM  code 
[33,  34,  35].  This  code  is  also  shock  fitting,  and 
solves  the  steady,  two  or  three  dimensional  Euler 
and  conservation  of  energy  equations.  A  cylindrical 
based  coordinate  system  transformed  to  a  compu¬ 
tational  space  is  used  to  allow  computations  on  an 
equally  spaced  grid.  The  difference  equations  are 
integrated  by  the  Beam- Warming  modified  upwind 
MacCormack  scheme  [36]  using  the  Van  Leer  flux 
vector  splitting  scheme  [37].  Characteristic  compat¬ 
ibility  relations  are  also  used  to  compute  the  surface 
pressure. 
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ORGANIZATION 


The  shock  shape  and  surface-pressure  distribu¬ 
tion  were  provided  to  the  Harris  finite-difference 
code  [38].  This  code  is  second-order  accurate,  and 
can  include  a  first-order  correction  for  variable  en¬ 
tropy.  In  the  first  iteration,  the  boundary  layer  is 
computed  as  if  the  stagnation  streamline  wetted  the 
entire  body.  Using  conservation  of  mass,  the  code 
then  traces  streamlines  upstream  from  where  they 
enter  the  boundary  layer  to  where  they  pass  through 
the  shock.  The  shock  angle  at  this  point  is  used 
to  compute  the  entropy  and  total  pressure  on  the 
streamline,  downstream.  In  the  second  iteration, 
this  variable  entropy  is  used  in  a  new  edge  condition, 
to  recompute  the  boundary  layer.  Although  it  is  well 
known  that  this  will  not  give  the  correct  boundary- 
layer  profile  in  the  outer  region  (Julius  Harris,  pri¬ 
vate  communication,  Sept.  1991),  this  code  is  fairly 
easy  to  use  and  readily  available.  This  code  runs  in 
a  few  seconds  on  a  400MHz  Pentium  II. 

Unfortunately,  the  computations  based  on  the 
first-order  Harris  code  showed  poor  agreement  with 
the  Stetson  blunt-cone  experiments,  as  will  be  dis¬ 
cussed.  Further  computations  are  to  be  carried  out 
using  a  viscous-shock-layer  (VSL)  code,  which  is  also 
readily  available  [39,  40].  This  method  solves  the  in- 
viscid  and  viscous  flows  in  a  coupled  way  that  is  sub¬ 
stantially  more  accurate,  while  remaining  substan¬ 
tially  simpler  and  faster  than  existing  Navier-Stokes 
or  Parabolized  Navier-Stokes  codes.  However,  VSL 
results  have  not  yet  been  obtained,  due  to  the  time 
required  to  learn  the  use  of  the  code.  This  VSL  code 
also  runs  in  a  few  seconds  on  a  400MHz  Pentium  II. 

The  boundary-layer  profiles  from  the  Harris 
code  were  then  supplied  to  the  eMALIK  code  for 
computation  of  instability  waves  [41,  42].  This  code 
solves  for  compressible  parallel-flow  linear  instabil¬ 
ity,  including  transverse  curvature  effects.  An  auto¬ 
mated  system  that  links  these  codes  was  previously 
developed  (e.g.,  Ref.  [43]).  Boundary-layer  profiles 
from  the  VSL  code  have  not  yet  been  used  for  sta¬ 
bility  analysis.  The  stability  code  can  run  several 
frequencies  along  the  length  of  a  cone  in  a  few  hours 
on  the  back  processor  of  a  clual-Pentium-II  400MHz 
PC.  This  speed  is  sufficient  for  design  purposes,  if  a 
suitably  validated  code  can  be  developed. 

Sharp  Cone  at  Mach  8 

General  Information 

Stetson  et  al.  carried  out  calibrated  measure¬ 
ments  of  instability  wave  growth  on  several  conical 
models  in  AEDC  Tunnel  B  at  Mach  8,  during  the 


late  1970’s  and  early  1980’s  [44,  45].  Detailed  hot¬ 
wire  measurements  were  carried  out  in  this  expen¬ 
sive  production  tunnel  by  J.  Donaldson,  under  the 
direction  of  K.  Stetson  [46,  47].  Both  are  now  re¬ 
tired,  although  Stetson  is  still  available  for  discus¬ 
sions.  The  experiments  were  focused  on  the  hyper¬ 
sonic  second-mode  instability,  which  is  likely  to  be 
dominant  on  smooth  convex  nearly-symmetric  ge¬ 
ometries  with  small  crossflow.  Although  these  are 
conventional-tunnel  measurements  with  high  ambi¬ 
ent  noise,  the  good-quality  detailed  measurements 
appear  to  deserve  further  computational  compar¬ 
isons.  Detailed  measurements  of  the  disturbances 
in  this  tunnel  are  available  [48]. 

The  model  was  a  7-degree  half-angle  cone  with  a 
40-inch  (1.016-m)  length  and  a  sharp  0.0015-inch  (38 
micron)  nose  radius.  The  cone  angle  of  attack  was 
zero,  to  within  the  accuracy  with  which  this  could 
be  set  and  measured.  The  model  was  in  thermal 
equilibrium,  but  heat  conduction  within  the  model 
was  not  negligible. 

Most  of  the  measurements  were  carried  out  at  a 
freestream  Mach  number  of  7.95,  with  a  total  pres¬ 
sure  of  225  psia  (1.55  MPa)  and  a  total  temperature 
of  1310  R  (728  K).  Additional  data  were  recorded 
at  a  lower  total  pressure,  to  determine  the  neutral- 
instability  region;  while  this  data  could  be  used  for 
linear-amplification  analysis,  the  amount  of  data  is 
limited  (see  Fig.  24  of  Ref.  [45]).  Data  was  also 
acquired  at  a  higher  pressure,  to  record  data  into 
transition  onset.  Thus,  transition-onset  data  should 
be  available. 

Measurements  were  carried  out  with  surface 
pressure  taps  (24)  and  thermocouples  (32).  Mean- 
flow  profiles  were  obtained  using  pitot-pressure  and 
total  temperature  probes.  Instability-wave  spec¬ 
tra  were  obtained  using  calibrated  constant-current 
hot-wire  anemometry.  Hot-wire  measurements  were 
recorded  at  one-inch  (2.54-cm)  intervals  between  10 
and  37  inches  (25.4  and  94.0  cm),  using  a  single  hot¬ 
wire  with  a  single  calibration.  Since  the  measure¬ 
ments  were  obtained  only  along  a  single  ray,  and  only 
at  the  height  above  the  wall  where  the  broadband 
amplitude  is  a  maximum,  Stetson  has  expressed  con¬ 
cern  about  three-dimensional  effects  (pp.  14-15  of 
Ref.  [45]). 

Growth  of  first  and  second- mode  instability 
waves  was  observed  under  cold-flow  conditions.  The 
second- mode  instability  was  dominant.  Second- 
harmonic  amplification  was  observed  downstream, 
a  good  indication  of  nonlinear  effects.  A  compari¬ 
son  of  the  second-harmonic  amplification  rate  to  a 
nonlinear  computation  would  be  interesting. 
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Although  Stetson  [44]  published  a  summary  of 
these  experiments,  extensive  tabulated  data  for  the 
mean  flow  remained  unpublished,  along  with  exten¬ 
sive  tabulations  of  the  instability  spectra  for  both 
uncalibrated  and  calibrated  fluctuations.  These  ta¬ 
bles  were  provided  courtesy  of  Roger  Kimmel  and 
have  been  entered  into  electronic  form  and  posted  on 
the  WG-10  website.  Many  detailed  measurements 
are  here  reported  for  the  first  time. 

The  amplification  rates  reported  by  Stetson  et 
al.  [44]  were  obtained  by  curve-fitting  and  differenti¬ 
ating  the  experimental  spectral-amplitude  data.  Al¬ 
though  most  of  the  past  comparisons  to  computation 
made  use  of  these  amplification  rates,  differentiation 
of  experimental  data  is  a  notoriously  uncertain  pro¬ 
cess  which  can  introduce  substantial  errors.  It  is  rec¬ 
ommended  that  future  comparisons  should  instead 
be  made  by  integrating  the  computed  amplification 
rates  and  comparing  amplitude  ratios  between  com¬ 
putation  and  experiment.  This  more-accurate  pro¬ 
cess  is  now  possible,  due  to  the  availability  of  the 
tabulated  spectral  amplitudes. 

Previous  Computational  Comparisons 

Mack  carried  out  the  first  set  of  detailed  com¬ 
parisons  to  this  experiment  [49,  50].  Three  meth¬ 
ods  were  used  to  obtain  the  mean-flow  profiles:  1) 
transformed  flat-plate  profiles,  2)  cone  profiles  from 
the  axisymmetric  boundary-layer  equations,  with¬ 
out  the  transverse  curvature  term,  and  3)  like  (2), 
but  with  the  transverse  curvature  term.  The  wall 
was  assumed  adiabatic.  Preliminary  computations 
including  the  effect  of  the  experimental  surface- 
temperature  gradient  are  reported  to  show  only 
minor  effects.  The  stability  was  analyzed  using 
parallel-flow  linear  instability  theory.  Locally  pla¬ 
nar  flow  is  assumed  in  the  stability  analysis.  Mack 
uses  the  Sutherland  viscosity  law  down  to  110. 4K, 
and  a  linear  law  below  this  (see  also  p.  57  of  Ref. 
[6],  where  this  law  is  stated  without  justification). 

The  effect  of  total  temperature  was  found  to  be 
significant.  The  effect  of  transverse  curvature  was 
small  (7%  for  one  amplification  rate).  At  R  =  1730, 
the  measured  frequency  of  the  maximum  amplifica¬ 
tion  rate  was  less  than  10%  higher  than  the  compu¬ 
tation.  However,  the  computed  amplification  rate  is 
about  50%  higher  than  the  measured  value.  Here, 
R  =  V ReXe  ■  the  square  root  of  the  usual  length 
Reynolds  number,  Rexe  =  Uex/ve.  Fig.  19  in  Ref. 
[44]  shows  that  the  experimental  amplification-rate 
curve  at  this  location  is  in  general  agreement  with 
earlier  measurements  by  Kendall  [51]  and  Demetri- 
ades  [52]. 


Fig.  5  in  Ref.  [50]  shows  that  the  computed 
maximum  amplification  rate  is  close  to  the  mea¬ 
sured  value  for  1200  <  R  <  1400,  but  the  measured 
value  increasingly  drops  below  the  computed  value 
for  R  >  1400.  Fig.  21  in  Ref.  [45]  shows  that  higher 
harmonics  also  become  significant  for  R  >  1400. 
Stetson  [45]  concludes  that  nonlinear  effects  become 
significant  for  R  >  1400,  making  comparisons  to  a 
linear  amplification  rate  inappropriate  there.  This 
is  a  very  important  point,  for  most  later  computa¬ 
tions  follow  Mack  in  comparing  to  the  experiment 
at  R  ~  1732.  Fig.  18  in  Ref.  [49]  shows  an  N-factor 
comparison  that  makes  the  same  point  -  the  exper¬ 
imental  amplification  rates  are  lower  than  the  com¬ 
puted  values,  particularly  above  R  ~  1800  where  the 
experimental  amplitudes  actually  decrease.  Mack 
suggests  that  this  is  due  to  nonlinear  effects  or  to 
the  start  of  transition. 

Mack  also  suggests  that  the  experimental  am¬ 
plification  rate  may  be  lower  because  a  pure  2D  nor¬ 
mal  mode  is  not  what  is  actually  being  measured. 
He  points  out  that  different  frequencies  may  have 
their  maximum  response  at  different  boundary-layer 
heights,  so  a  Fourier  analysis  of  the  signal  at  the 
height  where  the  broadband  amplitude  is  a  maxi¬ 
mum  does  not  necessarily  determine  the  single-mode 
amplification  rates. 

Chang  et  al.  [53]  show  one  figure  suggesting 
that  a  sharp-cone  PNS  solution  can  achieve  better 
agreement  at  R  —  1732,  but  few  details  are  given  in 
the  paper.  Their  linear-stability  equations  allow  for 
shock  motion.  A  similar  comparison  at  R  ~  1730 
is  shown  in  Fig.  4  of  Ref.  [54];  the  computational 
growth  rates  are  about  20%  above  the  experimental 
values.  Few  details  are  again  given. 

Chang  et  al.  [55]  performed  linear  and  nonlinear 
computations  using  PSE  methods.  The  mean  flow 
appears  to  be  obtained  using  the  Mangier  transfor¬ 
mation  from  the  flat  plate  solution,  but  the  details 
are  sketchy.  Chang  et  al.  also  find  poor  agreement 
between  linear-instability  theory  and  Stetson’s  am¬ 
plification  rates  at  R  ~  1730.  Nonparallel  effects  re¬ 
sult  in  an  increase  in  the  computed  growth  rates, 
making  the  agreement  worse.  However,  if  appropri¬ 
ate  nonlinear  effects  are  taken  into  account,  fairly 
good  agreement  can  be  obtained  (Figs.  9  and  13 
in  Ref.  [55]).  Chang  et  al.  suggest  that  nonlin¬ 
ear  effects  become  dominant  at  R  ~  1600.  Various 
assumptions  have  to  be  made  regarding  the  ampli¬ 
tude  of  the  particular  waves  that  are  included  in  the 
nonlinear  computation,  so  the  nonlinear  computa¬ 
tion  can  only  suggest  an  explanation  for  the  experi¬ 
ment,  it  cannot  be  validated  by  the  experiment.  The 
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linear-region  data  was  not  studied  in  detail,  proba¬ 
bly  because  Chang  et  al.  did  not  have  ready  access 
to  the  detailed  tabulated  data. 

A  series  of  detailed  computational  comparisons 
has  been  carried  out  by  Dallmann’s  group  at  DLR, 
Gottingen,  Germany.  Simen  et  al.  shows  com¬ 
parisons  carried  using  a  thin-layer  Navier-Stokes 
code,  combined  with  local  linear  instability  theory 
[56,  57,  58].  Fairly  good  agreement  is  obtained  for 
the  downstream  mean  profiles,  and  for  the  second¬ 
mode  growth  rates.  However,  these  comparisons 
were  later  repeated  by  Kufner,  using  the  same  code 
[59,  Section  6.2].  Kufner  did  not  obtain  the  same 
good  agreement,  even  when  using  Simen’s  grid.  Ac¬ 
cording  to  Stefan  Hein  (private  communications, 
October-December  2000)  and  Ref.  [60],  Kufner  be¬ 
lieved  that  Simen’s  good  agreement  was  due  to  insuf¬ 
ficient  convergence,  insufficient  grid  resolution,  in¬ 
sufficient  clustering  of  the  grid  near  the  shock,  and 
a  fortuitous  value  of  the  entropy-fix  parameter  that 
is  needed  to  stabilize  the  solution.  Kufner  obtained 
fully  converged  grid-independent  solutions  that  were 
independent  of  the  entropy- fix  parameter,  over  a  cer¬ 
tain  range.  These  solutions  resulted  in  amplifica¬ 
tion  rates  that  were  in  general  agreement  with  the 
other  computations,  and  that  disagree  with  the  ex¬ 
perimental  data.  Since  Esfahanian  also  obtained  for¬ 
tuitous  agreement  when  using  Simen’s  grid  (Stefan 
Hein,  private  communication,  October  2000),  it  ap¬ 
pears  that  Simen’s  good  agreement  was  a  misleading 
coincidence. 

In  summary,  at  the  present  time  it  appears  that 
Stetson’s  sharp-cone  experiments  show  waves  that 
rise  above  the  background  noise  for  R  ~  1200,  grow 
linearly  only  in  the  short  region  to  R  ~  1400,  and  be¬ 
come  nonlinear  downstream.  Most  detailed  compar¬ 
isons  to  linear  instability  have  mistakenly  been  car¬ 
ried  out  downstream,  where  the  growth  is  actually 
nonlinear,  and  comparisons  to  the  R  ~  1730  growth 
rates  therefore  show  poor  quantitative  agreement. 
Now  that  the  detailed  tabulated  data  are  available, 
new  comparisons  should  be  performed  for  the  linear 
instability  amplification. 


Mean  Flow  Comparisons 

Fig.  1  shows  the  wall  temperature  distributions 
measured  on  the  sharp  cone,  at  the  time  of  various 
mean-profile  measurement  runs.  The  coordinate  s  is 
the  arclength  along  the  surface.  The  x  locations  in 
the  legend  are  the  axial  locations  of  the  respective 
mean-profile  measurements,  and  Tw  is  the  surface 
wall  temperature.  The  forward  portion  of  the  model 
heats  up  with  time,  and  the  aft  portion  appears  to 
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Figure  1:  Wall  Temperatures  for  Stetson  Sharp  Cone 


cool  slightly.  The  data  do  not  agree  with  a  simple 
analysis  of  the  adiabatic  wall  temperature.  The 
temperatures  are  the  result  of  a  complex  heat  bal¬ 
ance  in  the  thin-wall  model,  which  incorporates  a 
solid  nosetip  and  a  water-cooled  sting  support. 

Fig.  2  and  3  show  the  tabulated  sharp  cone 
mean  profiles.  Both  show  overshoots  near  the 
edge  of  the  boundary  layer.  While  it  is  commonly 
accepted  that  the  total  temperature  can  overshoot 
near  the  boundary  layer  edge,  the  overshoot  in  the 
pitot  profile  is  suspicious. 

Fig.  4  shows  a  preliminary  comparison  to  a 
simple  computation.  The  comparison  is  shown 
for  Run  107,  at  Pt  =  225.3  psia,  Tt  =  1364. 4R, 
and  Mqo  =  7.95.  Here,  Pt  and  Tt  are  the  total 
pressure  and  temperature  in  the  stilling  chamber, 
and  Moo  is  the  freestream  Mach  number  in  the  test 
section.  The  measured  pitot  profiles  are  compared 
to  boundary  layer  computations  based  on  the  Har¬ 
ris  finite-difference  code  [38]  and  the  Taylor-Maccoll 
solution.  The  computational  profiles  were  reduced 
to  pitot  pressure  using  the  Rayleigh  pitot  formula 
above  Mach  1,  and  the  subsonic  isentropic  formula 
below  Mach  1.  The  streamwise  distance  x  =  35.0 
inches,  and  the  arc  length  distance  s  =  35.26  in. 
or  2.939  ft.  The  isothermal  computation  assumed 
a  1000R  wall  temperature,  which  is  near  the  mean 
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Figure  2:  Mean  Total  Temp.  Profiles  for  Stetson 
Sharp  Cone 


Figure  3:  Mean  Pitot  Profiles  for  Stetson  Sharp 
Cone 


Pt2,  meas.  pitot  pressure 


Figure  4:  Pitot  Profile  Comparison  for  Stetson 
Sharp  Cone 


of  the  approximate  measured  wall  temperature  dis¬ 
tribution.  The  adiabatic-wall  computation  results  in 
wall  temperatures  Tw  ~  0.856X),  while  the  measured 
is  Tw  ~  0.75 Tt. 

The  measurements  show  an  overshoot  in  the 
pitot  pressure  which  is  not  present  in  the  compu¬ 
tations.  Such  an  overshoot  has  been  observed  be¬ 
fore  for  blunt  cones  or  at  angle  of  attack  [61],  but 
the  author  is  not  aware  of  any  reliable  experimen¬ 
tal  data  showing  such  an  overshoot  for  a  sharp  cone 
at  zero  angle  of  attack.  The  obvious  experimental 
causes  might  include  small  residual  angles  of  attack 
or  a  bent  cone  tip.  It  was  also  thought  possible 
that  more  sophisticated  computational  approaches 
may  yield  this  overshoot.  However,  computations 
by  Kufner  using  the  thin-layer  Navier-Stokes  equa¬ 
tions  do  not  show  this  overshoot  in  the  pitot  pro¬ 
file  (Fig.  5).  These  computations  were  carried 
out  for  the  blunt-cone  conditions  of  Tt  =  750K, 
Tx  =  54.35K,  Moo  =  8.0,  Re x  =  8.202  x  106/m, 
7  =  1.4,  R  =  287.03J/kg-K,  and  with  a  Suther¬ 
land  constant  of  110K.  The  Mack  viscosity  law  is 
used  [6,  p.  57],  resulting  in  a  freestream  Reynolds 
number  that  is  6%  lower,  as  compared  to  the  use  of 
the  Sutherland  law  even  below  110K.  The  compu¬ 
tational  results  were  provided  by  Stefan  Hein  (pri¬ 
vate  communication,  April  and  October,  2000);  see 
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-e -  sharp  cone,  s=0.491 884  m 

blunt  cone,  s=0. 4871 93  m 


Figure  5:  TLNS  Pitot  Profiles  for  Stetson  Sharp  and 
Blunt  Cones  at  Blunt-Cone  Flow  Conditions,  from 
Kufner  via  Hein 

Ref.  [62].  Neither  the  blunt  cone  nor  the  sharp 
cone  (computed  at  the  blunt-cone  conditions,  not 
the  Stetson  sharp-cone  conditions)  shows  an  over¬ 
shoot  in  the  pitot  profile,  although  both  show  over¬ 
shoots  in  the  total-temperature  profile  (Fig.  6). 

However,  this  pitot  overshoot  is  similar  to  that 
observed  by  Kendall  on  a  flat  plate  when  insuffi¬ 
ciently  small  pitot  tubes  were  used.  Measurements 
with  smaller  pitot  tubes  did  not  show  an  overshoot. 
The  overshoot  was  then  attributed  by  Kendall  to 
probe  interference  [63].  A  similar  overshoot  is  ob¬ 
served  in  three  references  cited  in  Morkovin  et  al. 
[64];  Morkovin  also  explains  the  overshoot  as  probe 
interference.  Since  the  overshoot  decreases  with 
downstream  distance,  as  might  be  expected  since  the 
relative  pitot-probe  size  decreases  with  increasing 
boundary-layer  thickness,  it  seems  likely  that  probe 
interference  is  the  cause.  This  problem  then  casts 
doubt  on  all  the  mean-flow  profiles. 

Accurate  mean-flow  computations  at  the  Stet¬ 
son  sharp-cone  conditions  were  not  available  to  this 
author,  who  is  unaware  of  any  definitive  compar¬ 
isons.  However,  the  lack  of  a  pitot-profile  over¬ 
shoot  for  a  sharp  cone  was  confirmed  in  unpublished 
Navier-Stokes  computations  by  Manning  et  al.  (pri¬ 
vate  communication,  July  2000;  cp.  Ref.  [65]).  The 
pitot-profile  data  thus  appear  to  be  contaminated 


-e -  sharp  cone,  s=0.491 884  m 

-a -  blunt  cone,  s=0. 4871 93  m 


Figure  6:  TLNS  Total  Temperature  Profiles  for  Stet¬ 
son  Sharp  and  Blunt  Cones  at  Blunt-Cone  Flow 
Conditions,  from  Kufner  via  Hein 

with  probe  interference.  Since  the  total-temperature 
profile  data  were  obtained  at  the  same  time  with 
the  same  probe  assembly,  they  may  also  be  contam¬ 
inated;  however,  this  remains  to  be  shown  by  a  com¬ 
parison  to  accurate  computations. 

Inst  ability- Wave  Comparisons 

The  voltage  spectra  are  shown  in  Fig.  7.  Ev¬ 
ery  other  streamwise  station  is  shown,  to  improve 
clarity.  Since  the  measurements  were  made  at  the 
peak-amplitude  height  in  the  boundary  layer,  and 
since  the  sharp  cone  profiles  should  be  self-similar, 
all  measurements  should  be  at  the  same  mean-flow 
condition.  A  single  hot-wire  was  used  for  all  these 
measurements.  Thus,  all  local  calibrations  should 
be  identical,  and  ratios  of  voltage  amplitude  should 
be  ratios  of  flowfield  fluctuations.  However,  the  data 
remain  to  be  checked  in  detail.  The  variable  ZA  is 
the  hot-wire  height  above  the  wall,  in  inches,  and  s 
is  again  the  arc  length. 

The  calibrated  spectral  amplitudes  are  also 
available  and  are  published  here  for  the  first  time 
in  Figs.  8  and  9.  These  data  were  reduced  in  the 
early  1980’s  but  the  tabulated  data  were  not  pub¬ 
lished,  as  the  qualitative  features  do  not  differ  from 
the  uncalibrated  data.  The  printed  tables  were  re¬ 
duced  to  electronic  form  under  the  supervision  of 
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frequency,  Hz 

Figure  7:  Voltage  Spectra  for  Stetson  Sharp  Cone 


Figure  8:  Massflux  Spectra  for  Stetson  Sharp  Cone 


Roger  Kimmel  of  the  Air  Force  Research  Lab  (pri¬ 
vate  communication,  October  2000).  Figs.  8  and  9 
show  every  other  streamwise  station,  for  improved 
clarity.  The  run  numbers  and  measurement  heights 
are  the  same  as  in  the  uncalibrated  spectra.  Both 
spectra  are  normalized  by  the  mean  values  at  the 
probe  height.  Tabulated  data  are  available  from 
Roger  Kimmel  or  the  author,  or  from  the  WG  10 
website.  Uncalibrated  spectra  from  various 

freestream  hot-wire  measurements  are  also  available 
in  tabulated  form. 


The  integrated  growth  of  these  waves  was  com¬ 
puted  using  the  eMahk  code  and  compared  to  Mack’s 
results.  The  internal  similarity  solver  was  used  to 
obtain  profiles  for  the  eMallk  code.  The  Mack  data 
were  digitized  from  Fig.  6  in  Ref.  [49],  since  the  orig¬ 
inal  tabulated  data  are  lost  (Les  Mack,  private  com¬ 
munication,  July  2000).  Figure  10  shows  the  results, 
for  Tt  =  1310R,  Ree  =  1.42  x  106/ft.,  Me  =  6.8, 
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Figure  9:  Total- Temperature  Spectra  for  Stetson 
Sharp  Cone 
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Figure  10:  N-factor  Computations  for  Stetson  Sharp 
Cone;  Profiles  from  Internal  eMalik  Solver 


Pr  =  0.72,  and  an  adiabatic  wall.  Here,  Ree  and 
Me  are  the  unit  Reynolds  number  and  Mach  num¬ 
ber  at  the  boundary-layer  edge,  Pr  is  the  Prandtl 
number,  and  F  is  the  usual  nondimensional  fre¬ 
quency.  Significant  differences  exist;  in  particular, 
the  present  computations  find  no  instability  at  the 
higher  frequencies,  where  Mack  does  find  instability, 
and  the  N-factor  for  F  =  1.4  x  10-4  at  R  ~  2000 
differs  by  about  30%. 

The  cause  of  the  difference  is  difficult  to  deter¬ 
mine;  detailed  analysis  by  a  computational  expert  is 
clearly  required.  Possibilities  include  differences  in 
the  numerical  methods  and  the  fluid  property  mod¬ 
els,  such  as  the  heat-transfer  coefficient,  the  viscosity 
coefficient,  and  the  heat  capacity.  A  different  version 
of  the  eMahk  code  was  obtained  from  Roger  Kimmel 
(private  communication,  July  2000);  it  gave  results 
that  differed  substantially  at  high  frequencies  (e.g., 
from  N  =  1.0  to  N  =  1.4),  when  using  the  same 
input  file  and  when  executed  on  the  same  machine 
with  the  same  compiler.  The  present  version  was 
earlier  checked  against  Ref.  [42];  see  Ref.  [66]. 

A  set  of  computations  was  also  carried  out  using 
the  Taylor-Maccoll  solution  [67,  Sec.  16.5],  the  Har- 
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Figure  11:  N-factor  Computations  for  Stetson  Sharp  Figure  12:  N-factor  Computations  for  Stetson  Sharp 
Cone;  Profiles  from  Harris  Boundary-Layer  Code  Cone;  Four  Methods  for  One  Frequency 


ris  boundary-layer  code,  and  the  eMallk  code.  Figure 
11  shows  another  comparison  to  the  Mack  results, 
for  Tt  =  1310R,  Pt  =  225  psia,  Re oo  =  1.0  x  106/ ft., 
47oo  =  8.0,  Pr  =  0.72,  and  an  adiabatic  wall. 
Here,  the  Harris  code  was  modified  to  use  the  Mack- 
modified  form  of  the  Sutherland  viscosity  law.  The 
comparison  for  the  higher  frequencies  is  better  than 
that  obtained  using  the  internal  eMahk  solver,  due 
to  small  changes  in  the  boundary-layer  profiles  and 
their  derivatives  (not  shown) .  A  doubling  of  the  grid 
for  the  boundary-layer  solution  changed  the  present 
N-factor  results  by  less  than  0.5%. 

Figure  12  shows  a  final  comparison  between 
computations.  All  are  carried  out  at  Tt  =  1364. 4R, 
Pt  =  225.3  psia,  =  7.95,  and  for  an  adiabatic 
wall.  A  single  frequency  is  computed,  and  both  the 
N-factor  and  the  value  of  R  is  shown  vs.  the  stream- 
wise  distance  s.  For  the  first  two  sets,  the  profiles  are 
obtained  using  the  internal  eMallk  similarity  solver; 
for  the  last  two,  the  profiles  are  obtained  using  the 
Taylor-Maccoll  solution  and  the  Harris  boundary- 
layer  code.  The  first  pair  of  solutions  differ  in  the 
value  of  the  Prandtl  number  that  is  used.  The  first 
of  the  Harris-code  solutions  uses  the  Sutherland  vis¬ 


cosity  law,  even  below  110K;  the  second  uses  the 
Mack-modified  Sutherland  law.  For  the  first  3  so¬ 
lutions,  the  value  of  R  repeats  almost  exactly,  as 
it  should;  for  the  last  one,  R  differs  somewhat  due 
to  the  change  in  the  freestream  viscosity.  For  the 
two  solutions  using  the  internal  similarity  solver,  the 
small  change  in  Prandtl  number  changes  the  final  N 
factor  by  some  12%.  The  larger  of  these  two  solu¬ 
tions  is,  in  turn,  about  10%  below  the  Harris  code 
solutions.  The  change  in  viscosity  law  (for  the  Harris 
code  only)  changes  the  N  factor  only  slightly. 

The  best  method  of  comparing  to  the  exper¬ 
imental  data  is  to  integrate  the  computed  ampli¬ 
fication  factors  and  compare  the  amplitude  ratio 
to  the  experimental  value.  Since  the  experimental 
data  are  linear  only  to  about  R  =  1400,  and  since 
the  signal-to-noise  ratio  becomes  substantial  only  at 
about  R  =  1150,  the  comparison  was  carried  out 
from  R  ~  1150  to  R  ~  1350,  as  shown  in  Figure 
13.  The  experimental  voltage  spectra  were  used  for 
this  plot.  The  computation  was  carried  out  using 
the  Harris  boundary-layer  code,  for  Pt  =  225.3  psia, 
Tt  =  1364. 4R,  Mqq  =  7.95,  using  the  Sutherland  vis¬ 
cosity  law,  and  with  a  1000R  isothermal  wall  tem- 
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Figure  13:  Amplitude-Ratio  Comparisons  for  Stet¬ 
son  Sharp  Cone 


perature.  For  the  computation,  amplitude  ratios  are 
referenced  to  the  values  at  s  =  11.112  in.,  while  the 
reference  location  for  the  experiments  is  s  =  11.073 
in.,  which  is  sufficiently  close  for  this  initial  com¬ 
parison.  The  experimental  amplitude  growth  in  this 
region  is  small,  less  than  a  factor  of  3,  which  limits 
the  accuracy  of  any  comparisons.  The  early  non¬ 
linearity  is  probably  due  to  the  relatively  large  am¬ 
plitude  at  which  the  signal  becomes  larger  than  the 
noise  of  the  conventional  wind  tunnel.  Compare  Ref. 
[68],  where  growth  of  a  factor  of  2  is  measured  at  a 
far  lower  Reynolds  number  under  controlled  quiet 
conditions.  The  computed  maximum  amplitude  ra¬ 
tio  is  60%  larger  than  the  experimental  value,  and 
the  most  amplified  frequency  in  the  computation  is 
about  10%  lower.  The  experimental  peak  N  factor 
is  about  1.0,  while  the  computational  value  is  about 
1.5. 

Figure  14  shows  that  the  difference  is  probably 
not  due  to  the  effect  of  wall  temperature  distribu¬ 
tion.  Amplitude  ratios  are  shown,  as  in  Fig.  13, 
for  the  same  flow  conditions,  again  computed  using 
the  Harris  boundary-layer  code.  The  first  set  of  am¬ 
plitude  ratios  is  obtained  for  a  constant  isothermal 


Figure  14:  Amplitude-Ratio  Comparisons  for  Stet¬ 
son  Sharp  Cone:  Two  Wall  Temperature  Distribu¬ 
tions 

wall  temperature  of  1000R.  The  second  set  is  for  a 
linear  variation  of  wall  temperature  from  1060R  to 
966R;  this  approximates  the  experimental  distribu¬ 
tion  (Fig.  1).  The  reference  location  for  both  sets  of 
amplitude  ratios  is  s  =  11.112  inches.  The  effect  of 
the  varying  wall  temperature  is  small. 

Summary 

It  is  remarkable  that  accurate  agreement  be¬ 
tween  experiment  and  computation  remains  to  be 
obtained  for  the  growth  of  second-mode  waves  on  a 
round  sharp  cone  at  zero  angle  of  attack.  This  shows 
the  difficulty  of  accurately  measuring  and  comput¬ 
ing  hypersonic  instability.  The  experimental  mean- 
profile  data  is  clearly  contaminated  by  probe  inter¬ 
ference,  and  cannot  be  used  as  a  reliable  check  on 
the  mean-profile  computations.  This  means  that  it 
cannot  be  used  to  check  angle  of  attack,  mean  flow 
nonuniformity,  or  other  experimental  error  sources. 
The  early  nonlinearity  observed  in  the  experiments 
also  suggests  that  comparisons  of  linear  amplifica¬ 
tion  must  be  carried  out  either  with  controlled  dis¬ 
turbances  or  with  a  quiet  tunnel,  or  both. 

The  computational  results  shown  here  highlight 


1  - 12 


RTO-TR-AVT-007-V3 


Hypersonic  Laminar  Instability  on  Round  Cones  Near  Zero  Angle  of  Attack 


the  well-known  sensitivity  of  instability  computa¬ 
tions  to  small  errors  in  the  numerics  or  the  profiles. 
The  difficulty  no  longer  seems  to  be  computer  re¬ 
sources,  but  rather  a  robust  and  sufficiently  accurate 
code.  An  accurate  set  of  agreed-upon  results  should 
be  readily  available,  for  a  mean  flow  and  the  insta¬ 
bilities,  to  allow  benchmarking  to  be  carried  out.  It 
is  evident  that  this  is  needed  even  for  an  established 
code,  when  such  a  code  is  ported  to  a  new  location, 
user,  or  machine. 

The  poor  agreement  between  the  computation 
and  experiment  cannot  be  resolved  until  the  com¬ 
putational  and  experimental  methods  are  refined. 
The  present  experiments  also  suggest  that  nonlin¬ 
ear  effects  are  a  critical  part  of  second-mode  induced 
transition  on  sharp  cones,  since  so  much  of  the  wave 
growth  occurs  nonlinearly.  However,  this  remains  to 
be  confirmed  by  quiet-tunnel  experiments. 

Blunt  Cone  at  Mach  8: 

Stetson  et  al. 

General  Information 

In  addition  to  the  sharp  cone  measurements  dis¬ 
cussed  above,  Stetson  et  al.  carried  out  additional 
calibrated  measurements  of  instability  wave  growth 
on  several  blunt  cone  models  in  AEDC  Tunnel  B  at 
Mach  8  [69,  45].  The  apparatus  was  almost  identi¬ 
cal  [46,  47].  This  conventional  wind  tunnel  is  well 
characterized  [48].  The  blunt-cone  experiments  ex¬ 
amined  the  effects  of  the  entropy  layer  on  the  hyper¬ 
sonic  second-mode  instability. 

The  model  was  a  7-degree  half-angle  cone  with  a 
40-inch  (1.016-m)  length  and  a  0.15-inch  (3.81  mm) 
nose  radius.  The  cone  angle  of  attack  was  zero,  to 
within  the  accuracy  with  which  this  could  be  set 
and  measured.  The  model  was  in  thermal  equilib¬ 
rium  but  heat  conduction  within  the  model  was  not 
negligible.  Most  of  the  measurements  were  carried 
out  at  a  freestream  Mach  number  of  7.99,  with  a 
total  pressure  of  580  psia  (4.00  MPa)  and  a  total 
temperature  of  1350  R  (750  K). 

Measurements  were  obtained  with  surface  pres¬ 
sure  taps  (24)  and  thermocouples  (32).  Mean  pro¬ 
files  were  measured  using  pitot-pressure  and  to¬ 
tal  temperature  probes.  Instability  waves  were 
measured  using  calibrated  constant-current  hot-wire 
anemometry.  Hot-wire  spectra  were  recorded  at  one- 
inch  (2.54-cm)  intervals  between  10  and  31  inches 
(25.4  and  78.7  cm),  with  additional  measurements 
at  33,  35,  and  37  inches.  Four  hot-wires  were  used 
for  these  measurements,  each  calibrated  separately. 


There  is  no  data  for  transition  onset  on  the  blunt 
cone,  according  to  Stetson. 

Growth  of  the  first  and  second-mode  instabil¬ 
ity  waves  was  observed  on  a  symmetric  blunt  cone 
under  cold-flow  conditions.  The  second- mode  in¬ 
stability  is  dominant.  Although  amplification  was 
only  observed  downstream  of  the  nominal  swallow¬ 
ing  length,  strong  bluntness  effects  were  apparent. 
Some  entropy-layer  fluctuations  were  also  apparent, 
although  these  could  not  be  resolved  clearly  and 
no  quantitative  comparisons  have  been  attempted. 
Second-harmonic  amplification  of  the  second-mode 
waves  was  observed  downstream,  a  good  indication 
of  nonlinear  effects.  A  comparison  of  the  second¬ 
mode  amplification  rate  to  a  nonlinear  computation 
would  be  interesting. 

Previous  Computational  Comparisons 

Stetson’s  blunt-cone  flows  have  been  computed 
by  numerous  workers.  The  first  paper  that  is  known 
to  the  author  is  Malik  et  al.  [70],  who  used  a  Navier- 
Stokes  code  in  the  nose  region,  combined  with  a  PNS 
code  downstream,  followed  by  linear  instability  the¬ 
ory.  Most  computations  were  carried  out  for  the 
0.15-inch  nose  radius  that  is  described  here.  Ma¬ 
lik  et  al.  computed  a  peak  amplification  rate  at 
s/Rn  =  175  that  was  about  50%  higher  than  the 
experimental  value,  although  the  frequency  of  peak 
amplification  varies  by  only  a  few  percent.  Here,  Rn 
is  the  nose  radius.  They  speculate  that  the  differ¬ 
ence  may  be  due  to  nonlinear  effects,  since  a  second 
harmonic  is  present  there.  A  second  comparison  is 
shown  further  downstream  at  s/Rn  =  215;  the  peak 
amplification  rate  agrees  better  but  the  shape  of  the 
curve  is  substantially  different.  A  brief  description 
of  similar  results  is  also  shown  in  Ref.  [54]. 

Herbert  and  Esfahanian  et  al.  [71,  72]  carried 
out  a  detailed  comparison  to  the  blunt-cone  case. 
Fig.  10  in  Ref.  [72]  shows  that  the  actual  measured 
wall  temperatures  on  the  cone  are  about  20%  below 
adiabatic-wall  values.  Figs.  11  and  12  in  Ref.  [72] 
point  out  the  lack  of  repeatability  in  the  Stetson 
profile  data,  and  a  lack  of  agreement  in  the  profile 
shapes.  The  critical  importance  of  an  accurate  mean 
flow  computation  is  emphasized.  Fig.  15  in  Ref. 
[71]  shows  computed  values  of  the  peak  amplification 
rate  at  s/Rn  =  175  that  are  about  70%  above  the 
experimental  values,  although  the  frequency  of  peak 
amplification  agrees  to  within  a  few  percent. 

Kufner  et  al.  [60,  12]  carried  out  additional 
computations.  Figure  2  in  Ref.  [60]  summarizes 
the  comparisons  of  all  the  authors;  all  the  compu¬ 
tations  show  a  growth  rate  about  1.5  times  the  ex- 
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perimental  value,  at  a  frequency  roughly  10%  below 
the  experimental  value.  However,  these  comparisons 
were  carried  out  at  a  station  175  nose  radii  down¬ 
stream  of  the  tip.  Figure  7b  in  Stetson’s  1984  paper 
shows  substantial  amplification  of  second-harmonic 
frequencies  at  this  station.  Stetson  believes  that 
nonlinearity  is  substantial  at  this  station,  therefore 
a  growth-rate  comparison  to  linear  theory  has  dubi¬ 
ous  validity.  Many  workers  have  selected  this  station 
for  comparison  to  linear  theory,  following  Malik  et 
al.  [70],  and  apparently  through  insufficient  com¬ 
munication  with  Stetson  et  al.  Fig.  4  in  Ref.  [12] 
shows  a  significant  effect  of  the  model  for  the  second 
viscosity  of  air;  however,  this  effect  does  not  come 
near  to  accounting  for  the  discrepancy  with  experi¬ 
ment.  Fig.  31  in  Ref.  [12]  shows  that  computations 
based  on  a  second-order  boundary-layer  method  can 
give  good  agreement  with  computations  based  on 
the  thin-layer  Navier-Stokes  equations.  Ref.  [62]  re¬ 
ports  later  computations  at  the  Stetson  blunt-cone 
conditions,  but  does  not  provide  additional  compar¬ 
isons.  Ref.  [73]  reports  preliminary  comparisons  to 
the  entropy-layer  instability  reported  by  Stetson  et 
al.,  but  not  to  the  first  or  second-mode  instabilities. 

Mean  Flow  Comparisons 

Surface  temperature  distributions  are  shown  in 
Fig.  15.  These  were  collected  at  the  time  of  various 
mean- flow  profile  runs.  The  x  locations  in  the  leg¬ 
end  are  those  of  the  mean-flow  profiles.  In  the  Run 
101-105  series,  the  model  appears  to  heat  up  with 
time.  In  the  Run  74-76  series,  smaller  variations  are 
present.  Neither  dataset  agrees  with  a  simple  anal¬ 
ysis  of  the  adiabatic  wall  temperature. 

The  mean  flow  was  measured  for  this  exper¬ 
iment,  but  existing  comparisons  have  uncovered 
some  difficulties.  Fig.  5  in  Ref.  [71]  shows  that 
an  adiabatic-wall  computation  overpredicts  the  wall 
temperature  by  roughly  20%.  This  is  thought  to  be 
due  to  heat  transfer  to  the  water-cooled  sting  sup¬ 
port.  While  the  measurements  were  carried  out  af¬ 
ter  thermal  equilibrium  was  achieved,  heat  transfer 
within  the  model  apparently  cannot  be  neglected. 

Boundary-layer  profiles  were  also  measured  us¬ 
ing  pitot  and  total-temperature  probes,  as  shown  in 
Figs.  16  and  17.  The  total  temperature  plots 
show  runs  74-76  to  have  uniformly  thinner  profiles 
than  runs  101-104,  although  these  are  nominally  re¬ 
peat  runs. 

Unfortunately,  the  repeated  runs  reveal  discrep¬ 
ancies  of  10-15%  under  nominally  the  same  condi¬ 
tions  (see  also  Fig.  6  in  Ref.  [71]).  The  cause  of 
this  lack  of  repeatability  remains  uncertain.  Stetson 
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Figure  15:  Wall  Temperature  Distributions  for  the 
Stetson  Blunt  Cone 
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TTL,  total  temperature,  local,  deg.  Rankine 

Figure  16:  Mean  Total  Temp.  Profiles  for  Stetson 
Blunt  Cone 
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■a -  x=35,  run  1 01 


Figure  17:  Mean  Pitot  Profiles  for  Stetson  Blunt 
Cone 

was  always  concerned  that  these  profiles  may  also 
suffer  from  probe  blockage.  Boundary-layer  mea¬ 
surements  were  only  obtained  on  one  ray,  so  detec¬ 
tion  of  residual  three-dimensional  effects  is  limited 
to  surface  pressure  measurements  on  3  taps  at  90- 
deg.  azimuthal  intervals,  at  two  streamwise  stations. 
Although  these  generally  match  to  within  the  exper¬ 
imental  error,  pressure  is  an  insensitive  measure  of 
angle  of  attack. 

It  was  not  possible  to  repeat  many  runs,  due 
to  the  high  cost  of  operating  at  AEDC.  However, 
Stetson  believes  that  it  will  be  extremely  difficult  to 
obtain  higher  quality  data  in  a  production  facility. 
Furthermore,  budget  limitations  currently  prohibit 
further  tests  of  this  type. 

Fig.  18  compares  the  experimental  data  to  pro¬ 
files  computed  using  the  thin-layer  Navier-Stokes 
(TLNS)  equations  by  Kufner  et  al.  [12],  and  pro¬ 
vided  by  S.  Hein  (private  communication,  April 
2000).  The  conditions  are  as  described  with  Fig.  5, 
where  the  same  pitot  profile  is  shown;  an  adiabatic 
wall  is  assumed.  The  two  Stetson  datasets  are  re¬ 
peat  runs,  obtained  at  the  same  tunnel  conditions. 
The  Harris-code  computations  were  obtained  by  the 
present  author  for  Tt  =  1350R,  Pt  =  580  psia,  and 
Mx  =  8.00. 

Kufner’s  computation  shows  no  overshoot  in  the 
pitot  profile,  while  overshoots  are  observed  in  both 
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Figure  18:  Comparison  of  Pitot  Profiles  for  Stetson 
Blunt  Cone 

experimental  profiles.  Probe  interference  is  immedi¬ 
ately  suspected,  as  in  the  sharp-cone  case.  Although 
the  Kufner  computations  give  fair  agreement  for  the 
boundary-layer  thickness,  the  Harris-code  computa¬ 
tions  show  poor  agreement  for  both  the  edge  pitot 
pressure  and  the  boundary-layer  thickness. 

The  limitations  of  the  first-order  boundary-layer 
theory  are  thought  to  cause  the  problems  with  the 
Harris-code  results.  The  Harris  code  uses  a  variable 
VELEDG  to  define  the  edge  of  the  boundary-layer  as 
a  fraction  of  the  inviscid  edge  velocity.  This  param¬ 
eter  was  changed  to  99%  from  the  usual  value  of 
99.99%  (which  was  used  in  order  to  get  many  grid 
points  near  the  boundary-layer  edge  for  the  stability 
analysis).  Fig.  19  shows  the  results  for  the  same 
conditions  as  in  Fig.  18,  again  at  s  =  128.7Rn- 
The  edge  pitot  pressure  decreases  by  20%,  causing 
the  profiles  to  be  even  farther  away  from  the  ex¬ 
periment  and  the  more-accurate  Kufner  computa¬ 
tion.  This  effect  must  be  related  to  the  way  that  the 
first-order  boundary-layer  approach  obtains  the  edge 
conditions  by  tracing  streamlines  from  the  curved 
bow  shock.  The  use  of  a  first-order  boundary-layer 
method  for  computation  of  the  blunt-cone  profile  is 
clearly  inadequate  for  stability  purposes. 

Fig.  20  compares  total-temperature  profiles 


RTO-TR-AVT-007-V3 


1  -15 


Hypersonic  Laminar  Instability  on  Round  Cones  Near  Zero  Angle  of  Attack 


ORGANIZATION 


■a -  "edge"  =  99.99%  (usual) 


Figure  19:  Effect  of  Boundary-Layer  Edge  Parame¬ 
ter  on  Harris-Code  Profile 

from  the  same  Kufner  computation  to  the  Stetson 
profiles  from  two  repeat  runs.  The  Stetson  pro¬ 
files  were  obtained  at  s  =  127.7i?jv,  and  the  Kufner 
profile  is  from  s  =  127.9 R^.  Here,  the  three  pro¬ 
files  are  qualitatively  similar,  as  all  contain  simi¬ 
lar  overshoots  near  the  boundary-layer  edge.  How¬ 
ever,  the  quantitative  agreement  is  poor.  The  poor 
agreement  is  particularly  critical  when  one  remem¬ 
bers  that  the  stability  computations  are  highly  sen¬ 
sitive  to  the  second  derivatives  of  the  profiles.  Since 
Stetson’s  total-temperature  measurements  were  ob¬ 
tained  at  the  same  time  as  the  pitot  measurements, 
using  a  single  rather  large  rake,  both  pitot  and  total- 
temperature  measurements  may  be  contaminated  by 
probe  interference.  Also,  the  repeatability  of  the 
measurements  is  only  about  15%.  While  the  TLNS 
analysis  may  be  accurate  enough  for  stability  pur¬ 
poses,  the  quality  of  the  experimental  data  is  insuf¬ 
ficient  to  validate  the  computation. 

Inst  ability- Wave  Comparisons 

AEDC  carried  out  spectral  analyses  of  the  hot¬ 
wire  data  and  produced  tables  of  amplitude  vs.  fre¬ 
quency  at  these  streamwise  stations.  Tabulated  data 
is  available  both  for  the  uncalibrated  voltage  spectra 
and  for  calibrated  spectra  reported  for  total  temper¬ 
ature  and  massflux.  These  printed  tables  were  made 
available  by  Roger  Kimmel,  and  have  been  entered 


-  Kufner,  TLNS,  adiabatic  wall 

;ft  .  Stetson,  run  76 

■a -  Stetson,  run  104 


Figure  20:  Comparison  of  Total-Temperature  Pro¬ 
files  for  Stetson  Blunt  Cone 

into  electronic  form.  Figure  21  shows  the  massflow 
spectra.  Figure  22  shows  the  total  temperature 
spectra.  The  spectra  are  again  normalized  by  the 
mean  value  at  the  hot-wire  location,  given  by  the 
streamwise  arclength  s  and  the  probe  height  ZA.  The 
signal-to-noise  ratio  is  larger  for  the  total  tempera¬ 
ture  spectra,  since  the  constant-current  hot  wire  was 
operated  at  low  overheat. 

Rosenboom  et  al.  computed  second-mode  N- 
factors  for  this  case  [62],  using  an  accurate  Thin- 
Layer  Navier-Stokes  mean  flow.  The  instabilities 
were  computed  using  linear  local  and  nonlocal  meth¬ 
ods.  Hein  provided  the  tabulated  local-method  N 
factors  from  Fig.  17  of  Ref.  [62]  (private  com¬ 
munication,  August  2000).  These  data  are  replot¬ 
ted  in  Fig.  23,  and  compared  to  second-mode  N 
factors  computed  by  the  present  author  using  the 
Harris-code  boundary-layer  profiles  and  the  eMahk 
code.  While  the  110  and  120  kHz  results  are 
not  too  far  apart,  the  N  factors  for  the  higher  fre¬ 
quencies  differ  by  a  factor  of  more  than  2.  Both 
computations  were  carried  out  for  Tt  =  1350R  or 
750K,  Moo  =  8.00,  and  Rn  =  0.15  in.  or  3.81  mm. 
The  present  computation  uses  Pt  =  580  psia  and 
the  Sutherland  viscosity  law,  while  Rosenboom  et 
al.  used  Re oo  =  8.202  x  106/m,  which  is  obtained 
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Figure  21:  Massflux  Spectra  for  Stetson  Blunt  Cone 


Figure  22:  Total  Temperature  Spectra  for  Stetson 
Blunt  Cone 
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Figure  23:  N  factors  for  Blunt  Cone  Using  Two 
Methods 


from  this  same  total  pressure  by  using  the  Mack- 
modified  form  of  the  Sutherland  law.  If  the  nor¬ 
mal  Sutherland  law  is  used,  Re oo  =  8.76  x  106/m, 
which  is  only  7%  higher  and  insufficient  to  explain 
the  discrepancy.  The  substantial  errors  in  the  mean- 
flow  profiles  caused  by  the  use  of  the  first-order 
boundary-layer  code  result  in  very  large  errors  in  the 
highly-sensitive  instability  computation.  First-order 
boundary-layer  methods  are  inadequate  for  compu¬ 
tation  of  blunt-cone  instabilities,  even  for  prelimi¬ 
nary  design  purposes. 

Finally,  the  Kufner  N-factors  shown  in  Fig.  23 
were  used  to  compute  amplification  ratios,  and  com¬ 
pared  to  the  Stetson  amplification  ratios  for  total 
temperature.  The  total  temperature  data  was  se¬ 
lected  since  the  signal-to-noise  ratio  appears  bet¬ 
ter.  Run  56  at  s  =  32.254  in.  was  chosen  as  the 
terminal  station  for  a  linear  comparison,  since  Fig. 
22  shows  that  the  (nonlinear)  second  harmonic  be¬ 
comes  significant  further  downstream.  Run  59  at 
s  =  29.231  in.  was  chosen  as  the  reference  location, 
since  this  is  the  furthest  upstream  location  for  which 
the  amplified  frequencies  rise  clearly  above  the  back¬ 
ground  noise.  This  again  delineates  a  small  linear- 
amplification  region,  as  for  the  sharp  cone,  probably 
again  due  to  signal-to-noise  issues.  The  amplifica¬ 
tion  ratios  from  Run  56  to  Runs  58  and  59  were 
computed,  and  compared  to  values  that  were  lin¬ 
early  interpolated  from  the  Kufner  data  of  Fig.  23. 
Fig.  24  shows  the  results.  Both  methods  show  a 
small  amplification  of  about  1.4  or  1.5  to  S'  =  30.239 
in.,  and  both  show  a  larger  amplification  of  2.9  to 
3.3  to  S  =  32.254.  The  most  amplified  frequen¬ 
cies  differ  by  perhaps  10%  and  the  largest  amplifi¬ 
cation  factors  differ  by  perhaps  10%  also.  Outside 
the  range  of  computational  amplification,  the  exper¬ 
imental  amplification  factors  are  generally  not  reli¬ 
able,  due  to  signal-to-noise  problems.  While  the  fig¬ 
ure  again  shows  general  agreement,  it  is  still  clearly 
insufficient  for  reliable  code  validation.  Better  ex¬ 
perimental  data  is  needed. 


Summary 

Profiles  produced  by  a  first-order  boundary- 
layer  method  are  clearly  insufficient  for  stabil¬ 
ity  analyses,  although  second-order  boundary-layer 
techniques  do  appear  to  be  sufficient.  Good  agree¬ 
ment  between  computation  and  experiment  remains 
to  be  obtained.  Since  an  accurate  comparison  of  the 
mean  flow  profiles  is  not  possible,  due  to  probe  in¬ 
terference,  new  experiments  will  almost  certainly  be 
needed  to  resolve  the  disagreement. 
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Figure  24:  Amplification  Ratios  for  Blunt  Cone 

Sharp  Flared  Cone  at  Mach  6 

Three  instability  experiments  were  carried  out  in  the 
NASA  Langley  Mach-6  quiet  tunnel  before  it  was 
decommissioned  [74].  Comparisons  between  exper¬ 
imental  and  computational  N-factors  are  shown  for 
all  three  experiments  in  Fig.  11  of  Ref.  [74];  good 
agreement  was  obtained.  Due  to  schedule  limita¬ 
tions,  only  one  of  the  three  was  able  to  obtain  cali¬ 
brated  mean-flow  profiles  [75,  76,  77].  All  measure¬ 
ments  were  made  using  natural  fluctuations,  without 
controlled  disturbance  generators.  All  experiments 
were  carried  out  with  constant-voltage  anemome- 
try  (CVA),  due  to  electromagnetic  noise  problems 
with  more  conventional  anemometers.  No  system¬ 
atic  studies  of  CVA  calibration  techniques  had  been 
carried  out  at  that  time,  so  the  calibrations  are  lim¬ 
ited.  The  experiments  were  carried  out  on  both 
sharp  and  blunt  cones,  but  extensive  data  is  avail¬ 
able  only  for  the  sharp  cone. 

The  model  was  a  50.8-cm  sharp  cone.  The  first 
25.4  cm  formed  a  straight  cone  with  a  5-deg.  half¬ 
angle.  This  was  followed  by  a  gentle  flare  with  a 
2.364-m  radius  of  curvature.  The  profile  shape  is 
continuous;  only  the  second  derivative  is  discontin¬ 
uous  at  the  match  point.  The  tip  had  a  nominal 
radius  of  2.5  microns.  The  cone  angle  of  attack  was 


nearly  zero;  however,  most  measurements  were  car¬ 
ried  out  at  a  yaw  angle  of  approximately  0.1  deg. 
and  a  pitch  angle  of  approximately  0.1  deg.  [77,  p. 
20].  Appendix  D  of  Ref.  [77]  reports  additional  mea¬ 
surements  carried  out  after  additional  efforts  were 
made  to  align  the  model  at  zero  angle  of  attack. 

The  measurements  were  carried  out  at  Mach 
5.91,  with  a  total  pressure  of  896  kPa  and  a  total 
temperature  of  450K.  The  freestream  Reynolds  num¬ 
ber  was  9.25  x  106/m  (note  the  typographical  error 
onp.  2497  of  Ref.  [76]). 

Measurements  were  obtained  with  surface  pres¬ 
sure  taps  (29)  and  thermocouples  (51).  Mean  flow 
profiles  were  measured  using  hot  wires  and  CVA. 
Fluctuating  profiles  were  obtained  using  the  same 
hot  wires,  but  no  calibrated  fluctuations  were  re¬ 
ported.  Thus,  no  calibrated  eigenfunctions  are  avail¬ 
able. 

These  data  were  also  obtained  under 
equilibrium-wall  temperature  conditions.  That 
is,  the  measurements  were  obtained  only  after  the 
model-wall  thermocouples  reached  steady-state 
temperatures.  This  occurred  approximately  15-20 
minutes  after  initiation  of  Mach-6  flow,  following 
a  subsonic  preheat.  Unfortunately,  as  for  the 
Stetson  data,  this  is  not  the  same  as  adiabatic-wall 
conditions,  for  Chokani  states  that  the  heat-transfer 
within  the  model  is  not  negligible. 

Growth  of  first  and  second-mode  instability 
waves  was  observed  on  an  axisymmetric  cone  un¬ 
der  cold- flow  conditions,  in  a  quiet  tunnel  [78].  An 
adverse  pressure  gradient  exists  on  the  rear  half  of 
the  model.  Second- mode  amplification  is  observed. 
Gortler  interactions  are  possible  on  the  concave  flare, 
although  no  experimental  evidence  is  available. 

Lachowicz  et  al.  report  agreement  of  2-5% 
with  linear-instabilty  N-factor  computations  for  the 
second-mode  wave  growth  [76,  Fig.  9].  However, 
the  ratio  of  CVA  fluctuating  voltages  was  taken  as 
the  ratio  of  wave  amplitude.  This  assumption  was 
used  by  Stetson  and  others  for  constant-current  or 
constant-temperature  anemometry,  when  the  mea¬ 
surements  were  carried  out  with  the  hot  wire  at 
identical  mean-flow  conditions.  However,  the  valid¬ 
ity  of  this  assumption  for  the  Lachowicz  CVA  data 
is  presently  being  reviewed. 

If  this  assumption  can  be  verified,  and  any  errors 
caused  by  it  can  be  quantified,  this  dataset  may  be 
the  best  of  the  three.  Although  the  aft  end  of  the 
cone  was  outside  of  the  fully  quiet  region,  this  is  the 
only  data  obtained  under  quiet  conditions.  Static 
pressure  data  was  obtained  at  3  azimuthal  positions 
to  check  angle  of  attack,  although  all  boundary-layer 
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measurements  were  again  obtained  on  only  one  ray. 
The  mean  flow  profiles  are  in  fairly  good  agreement 
with  existing  computations. 

Additional  computations  have  now  been  carried 
out  (Ref.  [79]  and  [80]).  However,  the  experimental 
data  are  not  yet  available  in  tabulated  form;  detailed 
comparisons  may  be  carried  out  at  a  later  date.  New 
computations  were  reported  in  Ref.  [65],  but  wave 
amplifications  have  not  yet  been  compared. 

Summary 

Accurate  agreement  between  computation  and  ex¬ 
periment  remains  to  be  achieved  even  for  linear  in¬ 
stability  on  a  sharp  cone. 

The  boundary-layer  profiles  obtained  during  the 
Stetson  cone  experiments  suffer  from  probe  interfer¬ 
ence  and  are  also  repeatable  only  to  15%.  The  qual¬ 
ity  of  these  profiles  is  insufficient  for  code  validation. 
On  both  the  sharp  and  blunt  cones,  the  second-mode 
waves  amplify  by  less  than  a  factor  of  4  within  the 
linear  region,  after  achieving  a  significant  signal-to- 
noise  ratio.  This  is  insufficient  for  reliable  code  val¬ 
idation.  Conventional  tunnel  experiments  will  re¬ 
quire  controlled  perturbations  to  achieve  substan¬ 
tial  linear  amplifications  that  would  allow  an  accu¬ 
rate  validation  of  linear  computations.  It  is  expected 
that  a  linear  validation  will  usually  be  needed  before 
a  nonlinear  computation  can  then  be  validated. 

First-order  boundary-layer  methods  are  inad¬ 
equate  for  computation  of  blunt-cone  instabilities, 
even  for  preliminary  design  purposes. 
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